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Abstract 

A statistical mechanics argument relating partition functions to 
martingales is used to get a condition under which random geometric 
processes can describe interfaces in 2d statistical mechanics at critical¬ 
ity. Requiring multiple SLEs to satisfy this condition leads to some 
natural processes, which we study in this note. We give examples of 
such multiple SLEs and discuss how a choice of conformal block is 
related to geometric configuration of the interfaces and what is the 
physical meaning of mixed conformal blocks. We illustrate the general 
ideas on concrete computations, with applications to percolation and 
the Ising model. 
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1 Introduction 


Growth phenomena are ubiquitous around us. They have both very prac¬ 
tical applications and theoretical relevance. But they are rarely easy to 
study analytically and very few rigorous or exact results are known. In 
two dimensions, the description of a growing domain is often obtained in¬ 
directly through the description of a family of univalent holomorphic rep¬ 
resentations, leading quite generally to equations known under the name of 
Loewner chains. These techniques, based on the Riemann mapping theo¬ 
rem, are conceptually important but usually far from making the problem 
tractable. 

In the last few years, Loewner chains have been discovered which have 
a large hidden symmetry - conformal invariance - that makes them more 
amenable to an exact treatment These are known under the name of 
Stochastic or (Schramm) Loewner evolutions SLE. 

Their mathematical elegance and simplicity is not their sole virtue. They 
are also natural candidates to describe the continuum limit of an interface 
in two-dimensional statistical mechanics models at criticality. At the critical 
temperature and in the continuum limit, the system is believed to be con¬ 
formally invariant and physicists have developed many powerful techniques, 
known under the name conformal field theory or GET, to deal with local 
questions in a conformally invariant 2d system. However, nonlocal objects 
like interfaces posed new nontrivial problems that finally SLE could attack 
in a systematic way lilll. The connection between GET and SLE is now well 
understood iHiiainiiniiiaiii! and the interplay between the two approaches 
has proved fruitful. 

The way SLE describes an interface deserves some comments. As a 
guiding example, consider the Ising model in a simply connected domain, 
say on the hexagonal lattice. Suppose that the boundary is split in two arcs 
with endpoints say a and b and impose that on one arc the spins are up and 
on the other one the spins are down. In this situation each sample exhibits an 
interface. It joins the two points where the boundary conditions change and 
splits the domain in two pieces, one with all spins up on its boundary and 
one with all spins down. This interface fluctuates from sample to sample. 
What SLE teaches us is the following. Instead of describing the interface 
between a and b at once, SLE views it as a curve starting from say a and 
growing toward b. And SLE describes the distribution for the addition of an 
infinitesimal piece of interface when the beginning of the interface is already 
known. So the description is in terms of a growth process even if there was 
no growth process to start with. 
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As mentioned above, the probabilistic aspects of SLE as well as its con¬ 
nections with conformal field theory are now fairly well understood. How¬ 
ever, some fundamental questions remain, again directly related to natural 
questions in the statistical mechanics framework. 

The one we shall concentrate on in this note is what happens when, due 
to boundary conditions, the system contains several interfaces. Proposals 
for multiple SLEs have already been made in the literature nnnn, but our 
results point to a different picture. The simplest situation is in fact when 
there is only one interface but we want to deal with its two ends symmetri¬ 
cally so that two growth processes will interact with each other. Remember 
that standard SLE deals with the two ends of the interface asymmetrically. 
This has a price : time reversed SLE is an intricate object. 

As a guiding example for more than one interface, consider again the 
Ising model in a simply connected domain on the hexagonal lattice. If 
one changes boundary conditions from up to down to up and so on n = 2m 
times along the boundary, each sample will exhibit m interfaces, starting and 
ending on the boundary at points where the boundary conditions change, 
forming a so-called arch system. However, the interfaces will fluctuate from 
sample to sample and so does even the topology of the arch system. This 
topology, for instance, is an observable that is trivial for a single SLE. 

Our description will again be in terms of growth processes and Loewner 
chains. Eor standard SLE, the driving parameter is a continuous martingale 
and the tip of the curve separates two different states of the system (up and 
down spins for Ising), leading to a well defined boundary changing operator 
in statistical mechanics. The relation between the stochastic Loewner equa¬ 
tion and the boundary changing operator comes via a diffusion equation 
that they share in common. 

For multiple SLEs, we expect that for short time scales each curve grows 
under the inffuence of an independent martingale. At its tip stands the same 
boundary changing operator. But we also expect drift terms, describing 
interactions between the curves. 

The possibility of different arch topologies makes it even more natural to 
have a description with one curve growing at each boundary changing point 
so that each of them is on the same footing. So m interfaces are described by 
n = 2m growth processes of “half-interfaces” that finally pair in a consistent 
way to build arches. 

In statistical mechanics, each arch system has a well defined probability 
to show up. The law governing this finite probability space is again described 
by a Boltzmann weight which is nothing but a partial partition function. 

Our starting point is the reconsideration of the role of Boltzmann weights 
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and partition functions in statistical mechanics and their simple but cru¬ 
cial relationship with probabilistic martingales. This allows us to ask the 
question “by what kind of stochastic differential equations can one describe 
multiple SLEs ?” by imposing a martingale property and conformal invari¬ 
ance. This puts strong constraints on the drift terms and our main result 
is a description of the family of drift terms that are compatible with the 
basic rules of statistical mechanics. Each drift is expressed in terms of the 
partition function of the system. This partition function is given by a sum of 
Boltzmann weights for configurations that satisfy certain boundary condi¬ 
tions : at the starting points of the curves the boundary conditions change. 
The partition function depends on the position of these changes, so up to 
normalization, the partition function is in fact a correlation function. It 
satisfies a number of partial differential equations (one equation for each 
point) that are related to the diffusion equations for the multi SLE process. 
The solutions form a finite dimensional vector space. The positivity con¬ 
straint satisfied by physical partition functions singles out a cone which is 
expected, again guided by statistical mechanics, to have the same dimension 
of the underlying vector space and to be the convex hull of a family of half 
lines, so that a generic hyperplane section of the cone is a simplex. 

So geometrically, the drift terms are parametrized by a cone. Extremal 
drifts, i.e. drifts corresponding to extremal lines in this cone, lead to pro¬ 
cesses for which the final pattern formed by the growing curves is a given 
arch system. Drifts inside the simplex give rise to stochastic processes where 
the asymptotic arch system fluctuates from sample to sample. A crucial role 
to construct martingales describing interesting events is played by the short 
distance expansion in conformal field theory because this is what tells which 
terms in the partition function become dominant when an arch closes, i.e. 
when two driving processes of the multi SLE hit each other. 

The vector space of solution of the differential equations for the partition 
functions has a famous basis indexed by Dyck paths, which are in one to 
one correspondence with arch systems. But the basis elements do not in 
general correspond to extremal partition functions. We shall give a rationale 
for computing the matrix elements for the change of basis and compute a 
number of them, but we have no closed general formula. 

We shall illustrate our proposal with concrete computations for 1 to 3 
interfaces with applications to percolation and the Ising model. We shall 
also discuss the classical (deterministic) limit k ^ 0^, where only extremal 
drifts survive. 

The notes also cover the case when a number of boundary changes are 
very close to each other but the system is conditioned so that they do not 
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pair with each other. The details are in the main text. 

Our description is rather flexible in the sense that the speed of growth of 
each piece of interface can be tuned. Certain limiting cases lead to previously 
known processes which are examples of SLE(k, p) processes. 

It is appropriate here to stress that many of the probabilistic properties 
of the solutions of the stochastic differential equations that we introduce are 
conjectural at this point. We have made some consistency checks^ and the 
whole pattern is elegant, but our confidence comes more from our familiarity 
with conformal field theory and statistical mechanics. 

2 Basics of Schramm-Lowner evolutions: Chordal 
SLE 

Let us briefly recall what is meant by the chordal SLE — detailed studies 
can be found in [2j or |2j. The chordal SLE process in the upper half plane 
H is defined by the ordinary differential equation 

where the initial condition is gQ{z) = z and is a Brownian 

motion with variance parameter k > 0. Let < oo denote the explosion 
time of o with initial condition 2 ; and define the hull at time t hy Kt := 
{z G ]H[|r2 < t}. Then {Kt)t>o is a family of growing hulls, Kg C Kt for 
s < t. The complement H \ is simply connected and gt is the unique 
conformal mapping El \ iLt —> BI with gt{z) = z + o(l) at 2 ; —> 00 . One 
defines the SLE trace by 7 * = hm£|o 5 t~^(Ct + *e). The trace is a continuous 
path in El and it generates the hulls in the sense that EI\iLt is the unbounded 
component of El \ . Eor k < 4 the trace is a non-self-intersecting path 

and it doesn’t hit (9EI = M for t > 0 so LCj = 7 [o^t]. Eor 4 < /t < 8 a typical 
point 2 ; G El is swallowed, i.e. z & Kt for large t but 2 ; ^ 7[o,oo)- 
parameter range k > 8 the trace is space filling, 7 [o,oo) = 1®- Let us point 
out that no statistical mechanics models seem to correspond to k > 8. 

In the definition of chordal SLE we took the usual parametrization of 
time. Erom equation (Q) we see that gtiz) = z + 2tz~^ + 0{z~'^), which 
means (this could be taken as a definition) that the capacity of Kt from 

^For instance, Dubedat has derived general “commutation criteria” m for multiple 
SLEs. The processes we study are a special class satisfying commutation. This class 
extends vastly the special solution found by Dubedat, which in our language corresponds 
to self avoiding SLEs moving to infinity. 
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infinity is 2t. Since the capacity goes to infinity as t —> oo, the hulls Kt are 
not contained in any bounded subset of H. 

If the parametrization of time is left arbitrary, Schramm’s argument 
yields : 



where Ms is a continuous martingale with quadratic variation kQs (an in¬ 
creasing function going to infinity with s). In this formula, both qg and Mg 
are random objects. The capacity of Kg is 2qg. But this is not really more 


general than eq. © which is recovered by a random time change. 


3 A proposal for multiple SLEs 

The motivations for our proposal require a good amount of background, 
but the proposal and its main features themselves can be easily stated. We 
gather them in this section. Some of the results are conjectures. The rest 
of the paper will then be split into sections whose purpose will be either 
to motivate our proposal in general, or to prove its correctness in certain 
special but nontrivial cases by explicit computations. 

3.1 The basic equations 

We propose to describe the local growth of n interfaces in CFT, labeled 
by an integer i = 1, • • • ,re and joining fixed points on the boundary by a 
Loewner chain. We assume that 0 < k < 8 in the following. We list the set 
of necessary conditions and equations. 

Conformal invariance: The measure on nSLE is conformally invariant. 
Hence it is enough to give its definition when the domain D is the upper 
half plane H in the hydrodynamical normalization. 

Universe: The basic probabilistic objects are n (continuous, local) martin- 
gales Ml , z = 1, • • • , n with quadratic variation nqi ’ absolutely continuous 
with respect to dt and vanishing cross variation, defined on an appropri¬ 
ate probability space. By a time change we can and shall assume that 
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ferential equations 


dxP = dM« + Kdq? (a,, log Z) (Xf),... , Xf)) + ^ 


2dq- 


U) 




- X 


(i) 


( 2 ) 


(i) 

The initial conditions are Xq ^ = Xj ordered in such a way that Xi < X 2 < 

• • • < X„. 


Loewner chain: 

satisfies 


The map ft uniformizing the complement of the hulls 


dft{z) 


2dgf^ 

^ ft{z)-Xf' 


( 3 ) 


The initial condition is foiz) = z. With our conventions, the total capacity 
of the growing hulls at time t is 2t. 


Auxiliary function: The system depends on a function ,x„) 

which has to fulfill the following requirements : 

i) Z(xi, ■ ■ ■ , Xn) is defined and positive for xi < X 2 < • • • < Xn, 

ii) - ,Xn) is translation invariant and homogeneous. Its weight 
is hn- 2 m{i^) — nhi{K) for some nonnegative integer m < n/2, where^ 


2nhm{K) = m{2{m + 2) — k). 


in) Z{xi, ■ ■ ■ ,Xn) is annihilated by the n differential operators 




x-i — Xj 


-dr 




hl{K) 

{Xj - Xi)2 


We call this system of equations the nSLE system for n curves joining 
together the points Xi, - ■ ■ ,X„ and possibly the point at infinity. Systems 
for radial and dipolar versions of nSLE could be defined analogously. 


3.2 Arch probabilities 

It is known from GET that (relaxing the positivity constraint), the solutions 
to i), ii), Hi) form a vector space of dimension = (^) — (m-i) ~ 

(n + 1 — 2m) . 

The positive solutions form a cone and from the statistical mechanics 
interpretation, we conjecture that this cone has the same dimension and is 


more traditional notation for hm{K) is in the physics literature. 



generated by (i.e. is the convex hull of) dn,m half lines (extremal lines, pure 
states in the sense of statistical mechanics) so that a transverse section of 
the cone is a simplex. So each solution Z can be written in a unique way as 
a sum of extremal states. 

The numbers dn^m have many many combinatorial interpretations, but 
the one relevant for us is the following. Draw n + 1 points Xi < X 2 ■ ■ ■ < 
Xn < 00 ordered cyclically on the (extended) real line bounding the upper 
half plane H. Consider n — m disjoint curves in H such that each Xi is 
an end point of exactly 1 curve and 00 is an end point of exactly n — 2m 
curves. There are dn,m topologically inequivalent configurations, called arch 
configurations when n — 2m = 0. We keep the same name for m 7 ^ 0, writing 
archm configurations when precision is needed. 

Motivated by this, we claim the following ; 

a) To each arch configuration a corresponds an extremal state in 
the following sense : the solution of the nSLE system with Z oc Z^ can be 
defined up to a (possibly infinite) time, at which the growing curves have 
either paired together or joined the point at infinity and at that time the 
topology is that of the arch a with probability one. 

b) One can decompose a general solution Z of f), ii), in) as a sum of 

E 

a G axchm 


c) The probability that a solution of the nSLE system with auxiliary 
function Z ends in arch configuration a is the ratio 

Z^{Xi,--- ,Xn) 

Z{X,,--- ,Xn) 

evaluated at the initial condition (Xi, • • • ,Xn)- 

The hrst step toward a heuristic derivation of the above results will be to 
explain how to construct martingales - in particular martingales associated 
to interfaces - from statistical mechanics observables in a systematic way. 
But we start with a few comments. 
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4 First comments 


4.1 Statistical mechanics interpretation 

To have a specific example in mind, think again of the Ising model at the 
critical temperature. Let a be the lattice spacing. 

First, put n = 2m changes of boundary conditions from spins up to 
spins down and so on along the boundary at points xi/a,--- ,Xnla. In 
the continuum limit when a ^ 0 but xi,--- ,Xn have a finite limit, the 
partition function behaves like a homogeneous function Z{xi/a, ■ ■ ■ ,Xnla) 
of weight 0 (when both a and the Xi’s are rescaled) and CFT teaches us 
that Z{xi,--- ,Xn) satisfies i), ii), in) for n = 2m. Then, if a is an arch 
system, Z^ should be (proportional to the continuum limit of) the partial 
partition function when the sum of Boltzmann weights is performed only 
over the interface configurations with topology ol. 

To make generalized arch configurations, choose n and m with n > 2m. 
Put 2n — 2m changes of boundary conditions from spins up to spins down 
and so on along the boundary, n at points xi/o, • • • ,Xn/a and n — 2m at 
yi/a, ••• ,|/n- 2 m/a. Sum only over configurations where the interfaces do 
not joint two y-type points to each other. Take the continuum limit for the 
x’s as before, but impose that all y’s go to infinity and remain at a finite 
number of lattice spacings from each other. This is expected to lead again to 
a partition function Z(xi/a, • • • , Xnia) of weight 0 (when both a and the Xj’s 
are rescaled) and Z{xi, ■ ■ ■ ,Xn) satisfies i), ii), Hi) for the given n and m. 
If a is an arch^ configuration, Z^ should be (proportional to the continuum 
limit of) the partial partition function when the sum of Boltzmann weights 
is performed only over the interface configurations with topology a. 

Note that the prefactor between the continuum limit finite part and the 
real partition function is a power of the lattice spacing. The power depends 
on m, so it is likely to be unphysical to use a non homogeneous Z in the 
nSLE system, mixing different values of m for a fixed n. However, we shall 
later treat the example n = 2 mixing m = 0 and m = 1 because it is 
illustrative despite the fact that it breaks scale invariance. 

4.2 SLE as a special case of 2SLE 

For n = 2 the solution of i), ii), Hi) with m = 1 is elementary. Writing 
xi = a and X 2 = b, one finds Z (x. {b — Taking the first martingale 

to be a Brownian and the second one to be 0, one retrieves the equations for 
SLE growing from point a to point b in the hydro dynamical normalization. 
Let us recall briefly why. 


10 



We start from SLE from 0 to oo. The basic principle of conformal invari¬ 
ance makes the passage from this special case to the case when SLE goes 
from point a to point 5 on H a routine task. If u is any linear fractional 
transformation (i.e. any conformal transformation) from H to itself mapping 
0 to o and oo to b, the image of the SLE trace or hull from 0 to oo by u is 
by definition an SLE trace from a to 6 and this is measure preserving. The 
new uniformizing map is ht = u o o u~^ and it is readily checked that ^ 
is a rational function of ht whose precise form can be easily computed but 
does not concern us. 

Let us just mention that this rational function is regular everywhere 
(infinity included) except for a simple pole at ht = u{^t) and has a third 
order zero at ht = n(oo) = b. So the map ht is normalized in such a way 
that ht{b + e) = b + s + 0{e^), which is not the hydrodynamic normalization. 

But if Vt is any linear fractional transformation, Vt o gt o u~^ describes 
the same trace as ht = u o gt o u~^ . As long as the trace does not separate 
b from oo, i.e. as long as the trace has not hit the real axis in the segment 
]6, oo[, i.e. as long as oo is not in the hull, vt can be adjusted in such a 
way that ht = Vt o gt o u~^ is normalized hydrodynamically. Then dht/dt 
is a function of ht which is regular everywhere but for a single pole and 
vanishes at infinity, i.e. one can write dht/dt = 2/it/(/ij — at). The following 
computation is typical of the manipulations made with SLE (see e.g. @]). 
Write {gt o u~^){z) = w an compute from the definition 


Comparison gives 


dht , . du^ / ^ I V ^ ^ 

-^i^) = Vt{w) --. 

dt dt w — 


— {w) = 

dt ^ vt{w) - at w - it' 


But Vt is regular at w = it from which one infers that vt{it) = at (the poles 
in the two terms are at the same point) and /it = Vt(at)‘^ (the two residues 
add to 0). Going one step further in the expansion close to it yields ^(Ct) = 
-3v'/{it). Ito’s formula gives dat = -3v'/{it)dt + v't{it)dit + ^v'/{it)dt. So 
the time change /itdt = ds together with the definition dxs = Vti^t)dit yields 

dat(s) = dx^ + { k - 6) Ay^ds. 

But v'/(w)/v't{w)‘^ = 2/{vt{w) — vt{oo)) because vt is a linear fractional 
transformation. Einally, setting /it(s) = fs, ^t(s)(^) = ^t(s)(oo) = Bg and 
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vt{s)i^t{s)) = (^t{s) = -4^ we can summarize 


^ 2 

ds fs - As 


dBs 

ds 





dXs + (k - 6) 


ds 


- Bs 


where Xs is a Brownian motion with quadratic variation ns, fo = id, Aq = a, 
Bq = b. Thus chordal SLE from a to 6 in the hydrodynamical normalization 
is indeed a special case of 2SLE. 

The above equations are also a special case of SLE(fi;, p) {p = k — Q), but 
it should be clear that our general proposal goes in a different direction. 

As already mentioned, the description of chordal SLE from a to 6 in the 
hydrodynamical normalization in fact coincides with chordal SLE from a to 
b only up to the first time b is separated from oo by the trace. This time is 
infinite for k < 4, but it is finite with probability 1 for 4 < k < 8. The most 
obvious case is k = 6. The equation is nothing but the usual chordal SLEe 
ending at infinity, a consequence of locality (in the SLE sense, not in the 
quantum field theory sense used later). At that time, the real chordal SLE 
from o to 6 swallows oo, whereas the hydrodynamically normalized version 
swallows b. The solution to this problem is to use conformal invariance and 
restart the process again in the correct domain at the time when b and oo 
get separated by the trace. But this is not coded in the equations. 


4.3 Making sense 

The previous example should serve as a warning. Some serious mathematical 
work may have to be done even to make sense of our conjectures, let alone 
prove their correctness. The problems might be of different natures for k < 4 
and 4 < K < 8. We content with the following naive remarks. One of the 
problems is that the arches do not have to close at the same time. It may 
even happen that one of the growing curves touches the real line or another 
curve in such a way that the upper half plane is split in two domains and 
the one which is swallowed contains some of the growing curves. 

Our putative description of nSLE processes can be valid in this form only 
up the realization of such an event. The first thing to check should be that 
the event is realized with a probability obtained by summing Z^jZ over all 
a’s corresponding to compatible configurations (see figure^. In particular, 
the connected component of oo should contain at least noo > m — 1 curves 
for consistency, but that’s not an obvious property of our proposal. 

Consider the fate of the connected component of oo. If ni — m is even, 
conformal invariance suggests to continue the Loewner evolution simply by 
suppressing the points that have been swallowed, i.e. for the ni remaining 
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Figure 1: The probability of closing of an arch should be obtained by sum¬ 
ming ZajZ over all a’s corresponding to compatible configurations. Two 
compatible configurations are portrayed in the figure. 

points. If ni — m is odd, the same should be done, but the image of the point 
at which one interface has made a bridge should be included as a starting 
point for the continuation of the evolution. Preferably, the function Z for 
this new multiSLE system should not be adjusted by hand to make our 
conjectures correct, but should appear as a natural limit. We shall make 
comments on this and give concrete illustrations later. 

For the component that is swallowed, one can use conformal invariance 
again to change the normalization of the Loewner map in such a way that 
this component is the one that survives and then restart a new multiSEE for 
the appropriate number of points. This procedure may have to be iterated. 

Note also that our conjectures for arch probabilities do not involve any 

(i) 

details on the martingales . Indeed, we expect that there is some ro¬ 
bustness. But the precise criteria are beyond our understanding. 

4.4 A few martingales for nSLEs 

Our heuristic derivation of the nSLE system will in particular show that if 
Z also solves i), ii), Hi) (even relaxing positivity), the quotient 

ZjxP,--- 

is a local martingale. This can be proved directly using Ito’s formula. 

In particular, 

7 / ■v(l) 

is a local martingale bounded by 1, hence a martingale. On the other hand, 
a standard argument shows that if Pa is the probability that the system ends 
in a definite arch configuration a (once one has been able to make sense of 
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it) Pa{x[^\ ■ ■ ■ is a martingale. This is an encouraging sign. To get 

a full proof, one would need to analyze the behavior of Za{x[^\ ■ ■ ■ 
when one arch closes, or when one growing curve cuts the system in two, to 
get recursively a formula that looks heuristically like 

^a(^t 5 ‘ ‘ ‘ 5 ) r 

n 'l I'n'l ^ 

if the system forms asymptotically the arch system a' at large s. Such a 
formula rests on properties of Za(xi, ■ ■ ■ ,Xn) when some points come close 
together in a way reminiscent to the formation of arch a' : Z^'ixi, ■ ■ ■ ,Xn) 
should dominate all Z^s, a 7 ^ a' in such circumstances. In section|Hlwe shall 
use this to expand explicitly the Z^s in a basis of solutions to i), ii), in) 
which is familiar from CFT, very explicitly at least for small n. 


4.5 Classical limit 


Our proposal for nSLE has a non trivial classical limit at k —> 0+. The 

fO (i) 

martingales ' vanish in this limit, but the q) ^ remain arbitrary increasing 
functions. The function Z does not have a limit, but the Ui = ndxi log Z 
do. They are kind of Ricatti variables for which the equations read 


1 

2 



6 — 1 
2k, {xj — 


= 0 , 


which have a limit when k —> 0 ^, comparable to the classical limit of a 
Schroedinger equation. To summarize, the classical limit is 


dft{z) = 


2dq, 


(*) 


(d ■ 


dX« = U,iX^^\ 


Mz)-x) 

,xiYdqP+Y 


2dq, 


(J) 




xj"'> - X 


U)' 


where the auxiliary functions Ui{xi, ■ ■ ■ ,x„) are homogeneous functions of 
degree — 1 which satisfy dxi Uj = dx^ Ui and 



= 0 . 
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Figure 2: The factorisable Z leads to a very simple geometry. This case has 
been suggested previously with a slightly different approach. 

It is not too surprising that the differential equations for Z have become 
algebraic equations for the 17*’s, so that the space of solutions which was a 
connected manifold for k concentrates on a finite number of points in 
the classical limit. The classical system, maybe with an educated guess for 
the qi '’s, could be interesting for its own sake. 

4.6 Relations with other work 

Several processes involving several growing curves have appeared in the lit¬ 
erature. 

The first proposal was made by Cardy m- It can be formally obtained 
from ours by forgetting the conditions i), ii), Hi) and choosing a constant 
Z. The corresponding processes are interesting, but the relationship with 
interfaces in statistical mechanics and CFT is unclear for us. 

Dubedat m has derived a general criterion he calls commutativity to 
constrain the class of processes that could possibly be related to interfaces. 
Our proposal satisfies commutativity so they can be viewed as a special 
case satisfying other relevant physical constraints. Dubedat also came with 
a special solution of commutativity. It corresponds to the case m = 0 in 
our language. Then the space of solutions has dimension dnfl = 1 and the 
corresponding partition function is elementary: 

Z oc (4) 

i<j 

A single arch topology is possible, all interfaces converge to oo, see figure 
121 Maybe this is a good reason to call this case chordal nSLE. 
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5 CFT background 


There was never any doubt that SLEs are related to conformal field theories. 
The original approach [ZIIIIIIIIIIII used the operator formalism because if 
yields naturally martingale generating functions. Here, we use the correla¬ 
tor approach for a change. We restrict the presentation to a bare minimum, 
referring the newcomer to the many articles, reviews and books on the sub¬ 
ject f jl8[ll7| i. The reader who knows too little or too much about CFT can 
profitably skip this section. 

Observables in CFT can be classified according to their behavior under 
conformal maps. Local observables in quantum field theory are called fields. 
For instance, in the Ising model, on an arbitrary (discrete) domain, the 
average value of a product of spins on different (well separated) sites can be 
considered. Taking the continuum limit at the critical point, we expect that 
on arbitrary domains D there is a local observable, the spin. The product of 
two spins at nearest neighbor points corresponds to the energy operator. In 
the continuum limit, this will also lead to a local operator. In this limit, the 
lattice spacing has disappeared and one can expect a definite (but nontrivial) 
relationship between the energy operator and the product of two spin fields 
close to each other. As on the lattice the product of two spins at the same 
point is 1, we can expect that the identity observable also appears in such a 
product at short distances. Local fields come in two types, bulk fields whose 
argument runs over D and boundary fields whose argument runs over dD. 
In this paper, we shall not need bulk fields so we leave them aside. 

The simplest conformal transformations in the upper-half plane are real 
dilatations and boundary fields can be classified accordingly. It is customary 
to write (psix) to indicate that in a real dilatation by a factor A the field 
ips{x) picks a factor A*^. By a locality argument, boundary fields in a general 
domain D (not invariant under dilatations) can still be classified by the same 
quantum number. The number 5 is called the conformal weight of (ps- 

There are interesting situations in which (due to degeneracies) the action 
of dilatations cannot be diagonalized, leading to so called logarithmic CFT. 
While this more general setting is likely to be relevant for several aspects of 
SLE, we shall not need it in what follows. 

Under general conformal transformations, the simplest objects in CFT 
are so called primary fields. Their behavior is dictated by the simplest 
generalization of what happens under dilatations. Suppose psi r ' 'PSn 
boundary primary fields of weights , • • • ,Sn- If / is a conformal map from 
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domain L) to a domain D', CFT postulates that 

n n n 

(11 =(11 n 

j=i j=i j=i 

Symbolically, this can be written / : (ps{x) —> (ps{f{x))\f'{x)\^. It is interest¬ 
ing to make a comparison of these axioms with the previous computations 
relating chordal SLE from 0 to oo to chordal SEE from a to 6 in several 
normalizations. This also involved pure kinematics. 

As usual in quantum field theory, to a symmetry corresponds an observ¬ 
able implementing it. In GET, this leads to the stress tensor T{z) whose 
conservation equation reduces to holomorphicity. The fact that conformal 
transformations are pure kinematics translates into the fact that insertions of 
T in known correlation functions can be carried automatically, at least recur¬ 
sively. The behavior olT{z) under conformal transformations can be written 
as / : T{z) ^ T{f{z))f'{zf + c/l2Sf{z) where Sf = (/ 7 / 0 '- 1 / 2 (/ 7 / 0 ' 
is the Schwarzian derivative and c is a conformal anomaly, a number which 
is the most important numerical characteristic of a CFT. When c = 0, T is 
be a ( 2 , 0) primary field i.e. an holomorphic quadratic differential. When 
a (smooth) boundary is present, the Schwarz reflection principle allows to 
extend T by holomorphicity. Holomorphicity also implies that if O is any 
local (bulk or boundary) observable at point z £ D and v is vector field 
meromorphic close to z, the contour integral LyO = dwv{w)T(w)0 along 
an infinitesimal contour around z oriented counterclockwise is again a local 
field at z, corresponding to the infinitesimal variation of O under the map 
f{w) = w + sv{w). It is customary to write for v{w) = It is one 

of the postulates of CFT that all local fields can be obtained as descendants 
of primaries, i.e. by applying this construction recursively starting from 
primaries. The correlation functions of descendant fields are obtained in a 
routine way from correlations of the primaries. But descendant fields do not 
transform homogeneously. 

When V is holomorphic at x, LyO is a familiar object. For instance, if 
ifs is a primary boundary field, one checks readily that Lnfs = 0 for n > 1, 
Lq^Ps = dips and L^iips = [dxips]- The other descendants are in general 
more involved, but by definition the stress tensor T = L_ 2 /d is the simplest 
descendant of the identity Id. It does indeed not transform homogeneously. 

A primary field and its descendants form what is called a conformal 
family. Not all linear combinations of primaries and descendants need to 
be independent. The simplest example is the identity observable, which is 
primary with weight 0 and whose derivative along the boundary vanishes 
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identically^. By contour deformation, this leads to translation invariance of 
correlation functions when D has translation symmetry. 

The next example in order of complexity is of utmost importance for the 
rest of this paper. If {2h + l)c = 2h{8h — 5), the field 

—2{2h + + 2>L^_iLph 


is again a primary, i.e. it transforms homogeneously under conformal maps. 
In this case, consistent CFTs can be constructed for which it vanishes iden¬ 
tically. This puts further constraints on correlators. 

For example, when D is the upper half plane, so that the Schwarz prin¬ 
ciple extends T to the full plane, the contour for L _2 can be deformed and 
shrunken at infinity. Then, for an arbitrary boundary primary correlator 
one has 


2{2h + l) 


^1 + '^ 


a.=l 


1 _ Q __ 

Va-X (y« - xY 


{‘fs{oo) ipsc,iya)v’h{x)) = 0. (5) 


a=l 


It is customary to call this type of equation a null-vector equation. 

Note that the primary field of weight 5 sitting at oo has led to no con¬ 
tribution in this differential equation. Working the other way round, this 
equation valid for an arbitrary number of boundary primary fields with ar¬ 
bitrary weights characterizes the field (ph and the relation between h and 
the central charge c. 

The case of three points correlators is instructive. Global conformal 
invariance implies that 

{‘fs{y)‘fS'{y')‘Ph{x)) oc\y- - yf~^~^\y' - x\^~^'~'^. 


The proportionality constant might depend on the cyclic ordering of the 
three points. But if the differential equation for (ph is used, a further con¬ 
straint appears. The three point function can be non vanishing only if 

3(<5 - <5')^ - (2/i + l)(-5 + S') = h{h - 1). 

This computation has a dual interpretation ; consider a correlation function 
with any number of fields, among them a psiy) and a Phix)- If x and y 

^For other primary fields with the same weight if any, this does not have to be trne. 
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come very close to each other they can be replaced by an expansion in terms 
of local fields. This is called fusion. Several conformal families can appear 
in such an expansion, but within a conformal family, the most singular 
contribution is always from a primary. This argument applies even if c and 
h are arbitrary. But suppose they are related as above and the differential 
equation eq. © is valid. This equation is singular at x = y and at leading 
order the dominant balance leads to an equation where the other points are 
spectators. One finds that the only conformal families that can appear are 
the ones whose conformal weight 5' satisfies the fusion rule. 

This is enough CFT background for the rest of this paper. We are now 
in position to give the heuristic argument that leads to our main claims. 

6 Martingales from statistical mechanics 

The purpose of this section is to emphasize the intimate connection between 
the basic rules of statistical mechanics and martingales. The connection 
is somehow tautological, because statistical mechanics works with partition 
functions, i.e. unnormalized probability distributions, all the time. In the 
discrete setting, this makes conditional expectations a totally transparent 
operation that one performs without thinking and even without giving it a 
name. But the following argument is, despite its simplicity and its abstract 
nonsense flavor, the crucial observation that allows us to relate interfaces in 
conformally invariant statistical mechanics to SLEs. 

6.1 Tautological martingales 

Consider a model of statistical mechanics with a finite but arbitrarily large 
set of possible states S. Usually one starts with models defined on finite 
grid domains so < oo is natural. To each state s £ S we associate a 
Boltzmann weight® w{s). The partition function is Z = so that 

it normalizes the Boltzmann weights to probabilities, P{s} = Since 

S is finite, we can use the power set V{S) = {U : U C 5} as a sigma 
algebra. The expected value of a random variable O ; S' —> C is denoted by 
E[0] = {0) = ^ZsesO{s)w{s). 

Note that if {Sa)a£i is a collection of disjoint subsets of S such that 
Ua£iSa = S, then the collection of all unions T = {Uag//SQ : C /} is a 

^Usually the Boltzmann weight is related to the energy It(s) of the state s through 
to(s) = exp(—/3tif(s)), where d is the inverse temperature (a Lagrangian multiplier related 
to temperature, anyway). 
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sigma algebra on S. Conversely, since S is finite, any sigma algebra ^ on S' 
is of this type. 

Consider a filtration, that is an increasing family of sigma alge¬ 

bras {0, S} C C JCj c V{S) for all 0 < s < t. Denote the corresponding 
collections of disjoint sets by and define the partial partition func¬ 

tion Za^ = V (t) w{s). The conditional expectation values 


{0)t 


aGit 


0{s)w{s) 


1 


s 


(t) 


E (i E 


(t) 


are martingales by definition: for s < t we have 

E[E[0\J^t] \J^s] =E[0\J^,] 

Notice that the probability of the event S® is conveniently P[Sa^] = jZ. 

Suppose that the model is defined in a domain D G C and that there are 
interfaces in the model. Parametrize portions of these interfaces touching 
the boundary by an arbitrary “time” parameter t in such a way that n paths 
i = 1 , • • • ,n (which are pieces of the random interfaces) emerge from 
the boundary at t = 0 and are disjoint at least when t is small enough, 
see figure |31 To avoid confusion we write the time parameter t now as a 
subscript and continue to indicate the dependence of s e 5 by parenthesis, 
so 11 -^ 7 ^^*^ (s) is a parametrization of the piece of interface if the system 
is at state s. Then we can consider the natural filtration of the interface 
by taking = <^( 74 *^ : 0 < t' < t, i = 1 , • • • , n) to be the sigma algebra 
generated by the random variables 7 ), up to time t. 

The boundary conditions of the model are often such that conditioning 
on the is the same as considering the model in a smaller domain (a 

part of the interface removed) but with same type of boundary conditions. 
Of course the position at which the new interface should start is where the 
original interface would have continued, that is the 7 ) "^’s. Let Dt be the 
domain D with the 7 jq 7 removed. 

The starting point of the next section is the input of conformal invariance 
in this setup. 
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Figure 3: A discrete statistical mechanics model with portions of interfaces 
specified. 


6.2 Simplifying tautological martingales 


We start from the situation and notations at the end of the previous section. 
If in addition we are considering a model at its critical point, then the 
continuum limit may be described by a conformal field theory. At least for 
a wide class of natural observables O, the expectation values become CFT 
correlation functions in the domain D of the model 


(O) 


T.s&s^is)w{s) 

Z 


^^)CFT, b.c. 

/..vCFT, b.c. 
A/D 


We need to write the correlation function of identity (proportional to Z) in 
the denominator because the boundary conditions (b.c.) of the model may 
already have led to insertions of boundary changing operators that we have 
not mentioned explicitly. 

The closed martingales become 


(o>< = E ^ E 


E—/ y(t) 
a&h 


/^\CFT, b.c. 
W/Dt 
/..vCFT, b.c. 


where in the continuum limit Dt might be D with hulls (and not only traces) 
removed. 

For certain (but not all) observables, {O) is computing a probability, 
which in a conformal field theory ought to be conformally invariant. But 
{O) is written as a quotient, and this means that the numerator and denom¬ 
inator should transform homogeneously (and with the same factor) under 
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conformal transformations. In particular, the denominator should transform 
homogeneously. This means that (1)^^"*"’ - which depends on the posi¬ 

tion of the boundary condition changes - behaves like a product of boundary 
primary fields. Then, by locality, for any O, the transformation of the nu¬ 
merator under conformal maps will split in two pieces: one containing the 
transformations of O and the other one canceling with the factor in the 
denominator. So we infer the existence in the CFT of a primary boundary 
field, denoted by 'ijj{x) in what follows, which implements boundary condi¬ 
tion changes at which interfaces anchor. Hence we may write 


and 

^^)CFT, b.c. ^ 

As will become clear later, there might also be one further boundary 
operator anchoring several interfaces. We do not mention it explicitly here 
because it will sit at a point which will not be affected by the conformal 
transformations that we use. 

Write the transformation of the observable O as / : O —> under a 

conformal map. Denote by ft a conformal representation ft'.Dt^D and 
write /(7i*^) = The expression for the closed martingale {0)t can now 
be simplified further 




( 6 ) 


The Jacobians coming from the transformations of the boundary changing 
primary field if have canceled in the numerator and denominator. The 
explicit value of the conformal weight of if does not appear in this formula. 

Of course, we have cheated. For the actual map ft which is singular at 
(i) 

the 7 ^ '^s these Jacobians are infinite. A more proper “derivation” would go 
through a regularization but locality should ensure that the naive formula 
remains valid when the regularization is removed. Eq. m is the starting 
point of our analysis. 


7 Derivation of the proposal 

The heuristics we follow is to describe a growth process of interfaces by a 
Loewner chain ft compatible with conformal invariance in that the right 
hand side of eq. is a martingale. 
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7.1 The three ingredients 


Loewner chain: If we use the upper half plane as a domain, D = M, and 
impose the hydrodynamic normalization, the equation for ft has to be of 
the form 


dft{z) 


i /.(^) - 4” 


for some processes X 


(4 


i = 1,... ,n. The initial condition is fo{z) = z. 


Interfaces grow independently of each other on very short time 
scales: Schramm’s argnment deals with the case of a single point. We expect 
that on very short time scales the growth processes do not feel each other 
and Schramm’s argnment is still valid, so that dX^' = dM^ F^' where 

the Mj*^’s are n (continuous,local) martingales with quadratic variation Kqt^^ 

(i) 

and vanishing cross variation and T) is a drift term. 

The martingale property fixes the drift term: The drift term will be 
computed by imposing the martingale condition on {0)t when O is a product 
of an arbitrary number I of boundary primary fields O = 0^=1 
The insertion points are away from the bonndary changing operators and ft 
is regular with positive derivative there. Substitution of in formnla ® 
yields 


<11 »•- (nt,^(4b)g- ’ ■ 


a=l 


a=l 


(7) 


7.2 Computation of the Ito derivative of (fl 

In formula 0 , denote respectively by Zf, Zt and Jf the nnmerator, de¬ 
nominator and Jacobian factor on the right hand side. 

It is useful to break the computation of (y(")))j in several 

steps. 


- Preliminaries. 

Ito’s formnla for the '!/>’s gives 


dV>(X«) = V^'(X«)(dM« + F«) + ^r{X?)dqf\ 
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Using the Loewner chain for ft{z) and its derivative with respect to z, one 
checks that 


i 


( MtjY)) 

V/i(y)-x« 


6Mft{Y)) \ 


- The Ito derivative of Zfjf. 

(i) 

The time t being given, we can simplify the notation. Set xt = and 
Ua = /t(y^“^) and apply the chain rule to get 


dizfjr^ 


TV 





- First use of the null-veetor equation : identification of ip. 

Let us concentrate for a moment on the familiar chordal SLE case, for which 
n = 1. The drift term is known to be zero. The boundary conditions 
also change at oo (the endpoint of the interface) and there is an operator 
there, that we have not written explicitly because the notation is heavy 
enough. Anyway, Zt is a two-point function with one of the fields at infinity, 
so it is a constant. For chordal SLE, the drift term in the Ito derivative of 
the putative martingale vanishes if and only if 



+2E 


Va - X 


f) __ 


yf = 0, 


where for simplicity we have written x = xi. 

Comparison with eq. © implies that ip has a vanishing descendant at level 
two and has conformal weight h = = /ii(k) : 


i’ix) = Thi{K){x)- 

The central charge is c = . 

This is of course nothing but the correlation function formalism version of 
the original argument relating SLE to GET, which was given in the operator 
formalism, see 

- Second use of the null-veetor equation. 

Now that Ip has been identified, we can return to the general case, with an 
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arbitrary number n of growing curves. Each growing curve has its own field 
Tp and each field ip comes with its differential equation, which is eq.Q but 
for / + n — 1 spectator helds, the I fields (p and the n — 1 other ip's. So Zf 
is annihilated by the n differential operators 




1 _ Q __ 

{ya-x,)P 


_ya - Xi 

We can use this to get a simplihed formula 


+2E 


1 ^ /ii(k) 

^^3 ( \2 

2 ) 


d{zfjf) = Jf VZf , dZt = VZt 


where V is the first order differential operator 


E 






1 

Xj — Xi 


dfr 


{Xj-X,)^\ j 


The formula for Zt is just the special case I = 0. 
- Final application of Ito’s formula. 



where Q is the first order differential operator 


E 




KdqiP\d^, log Zt)-2 

j¥=i 




U) 


Xi — Xi 


d. 


Xi 


The martingale property is satisfied if and only if the drift terms vanish. 


7.3 Main claim 

To summarize, we have shown that the system 


= Y 


2dq. 


(i) 




id 


, dX® = + F^^ 


admits conditioned correlation functions from CFT as martingales if and 
only if 


F,^^ = Kdqf\di,^ logZt) + 2Y 


Pd 




U) 
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where Zt is a partition function. It is under this condition that it describes 
the growth of n interfaces in a way compatible with statistical mechanics 
and conformal field theory. 

In fact, we have used a special family of correlators. But the same 
argument applies to all operators (hence the “if” part). Of special interest 
in the sequel will be the case when O is a topological observable, for instance 
taking value 1 if the interface forms a given arch system and 0 otherwise. 
No Jacobian is involved for such observables and the numerator looks again 
like a partition function. 


7.4 The moduli space 

From the definition of Zt as a correlation of primary fields with null descen¬ 
dants at level 2, it is clear that properties i), ii), in) are satisfied, except 
maybe for the quantization of the possible scaling dimensions of Zt, to which 
we turn now. 

This is standard material from CFT and we include it here for complete¬ 
ness. 

The correlator (oo)V'(xi) • • • V’(a:n)) on the real line satisfies n dif¬ 
ferential equations. We shall recall why the space of simultaneous solutions 
which have global conformal invariance has dimension 


n 


m 


n 


= {n + 1 — 2m)- 


n\ 


^m — lj ' ' ^ ml{n — m + 1)\ 

if hoo = hn- 2 m{i^) for some nonnegative integer m < n/2 and has dimension 
0 otherwise. This will end the derivation of our proposal and match the 
counting of arches. 

At the end of the background on conformal invariance, we mentioned 
fusion rules; when ^hi{K){x) and a ^hj{K){y) are brought close together, 
they can be expanded in a basis of local operators that can be grouped in 
conformal families. We also recalled why the weight h' of the primaries in 
each conformal family had to satisfy 2>{hj{K)—h')‘^ — (2hi{K)-\-l){hj {K)-\-h') = 
/ii(k)(/ii(k) — I), so that only two conformal families can appear in a fusion 
with (pht(K)- The two conformal weights are easily found to be h' = hj±i{K). 
Furthermore, = 0 and one can show that the corresponding field has 

to be the boundary identity operator. By global conformal invariance, the 
only local operator with a non vanishing one point correlator is the identity 
and boundary two point functions vanish unless the two local fields have the 
same conformal weight. This takes care of the counting and selection rules 
for the n = 0,1 cases. 
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One proceeds by recursion. The points are ordered xi < X 2 • • • < 

If n > 2 then move X 2 close the xi (for instance by a global conformal 
transformation) and fuse to get an expansion for local fields at xi say. Only 
the conformal families of ^hi±-i_{K) appear. If n = 2 this fixes the weight of the 
field at oo. If n > 3, iterate. This leads immediately to the selection rules 
mentioned above : the field at infinity has to be a dimension 

is nothing but the number of path of n steps ±1 from 0 to n — 2m on 
the nonnegative integers, a standard combinatorial problem whose answer 
(m) “ (m-i) ■ efficient way to do the counting is by the reflection 

principle. The possible outcomes of each fusion can be encoded in a so- 
called Bratelli diagram: 


/i4(K) 


hl{K) 


/12(k) 


\ / 


/ io ( k ) ho{K) 

ISLE 2SLE 


hsin) 


hi{K) 


3SLE 


h2{K) 

hoin) 

4SLE 


( 8 ) 


This is totally parallel to the discussion of composition of n spins 1/2 for 
the representation theory of the Lie algebra of rotations. The multiplicity 
is exactly one when m = 0 which corresponds to the partition function @ 
and to the insertion of the operator ^hn{K) infinity, toward which the n 
interfaces run. 

What is not proved here is that the different paths lead to a basis of 
solutions of the n partial differential equations, but it is true. Each path 
corresponds to a succession of choices of a single conformal family, one at 
each fusion step. Let us mention in advance that multi SLE processes, i.e. 
the consideration of multiple interfaces, will lead to the definition of another 
basis with a topological interpretation. 
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8 Multiple SLEs describing several interfaces 


8.1 Double SLEs 

The case of double SLEs is instructive and simple to analyze. Although 
double SLEs is sometimes interesting for its own sake, the purpose of this 
section is to give easy examples to guide the study of the general case. 


8.1.1 2SLEs and Bessel processes 

To specify the process we have to specify the partition function Z. There are 
only two possible choices corresponding to two different type of boundary 
conditions, or alternatively to two different fields inserted at infinity; 

(hoo|V’(^i)V’(^2)|0) = const. X (Ai - A 2 )'^ 

where the exponent is A = Loo ~ 2/ii(«;) and the constant will be fixed 
to 1 from now on. According to GET fusion rules. Loo can only be either 
^ 2 (k) = or Lo(k) = 0. The exponent becomes A = 2/«;orA = 
respectively, so that we have two basic choices for Z; 

Zo = (Ai - or Z 2 = iXi-X 2 f/^ 


As we shall see, choosing Zq selects configurations with no curve ending at 
infinity - so that we are actually describing standard chordal SLE joining 
to the two initial positions of Xi and X 2 - while choosing Z 2 selects config¬ 
urations with two curves emerging from the initial positions of Xi and A 2 
and ending both at infinity. 

Up to normalizing the quadratic variation by dql = aidt so that the 

( 2 ) ( 2 ) ( 2 ) 

martingales are simply dM^ ^ = yJnaidBi ^ with dBl two independent 
normalized Brownian motions, our double SLE equations become : 


dftiz) ^ 


dXi'l = 




2ai dt 


+ 


202 


,- ,„(!) , 2a2 + KAai 

V^dB, - ^dt 

Xt - x^. 

^ dBf^ + 


- X 


( 1 ) 


It describes two curves emerging from points Xi = Xg^^ and X 2 = X^'^ at 
speeds parametrized by oi and 02 . 
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Figure 4: The two geometries for 2SLE: on the left is the ease hoo 
on the right hoo = h 2 {K). 


0 and 


Up to an irrelevant translation, the process is actually driven by the 
difference Yt = Up to a time change, ds = K(ai + a 2 )dt, this is 

a Bessel process, 

dYs = dBs + ^ ds, 

of effective dimension deff = 1 + 2A + 4/k. For hoo = ^ 2 (k) (he. A = 2/k) 
the dimension is deff = 1 + 8 /k and for hoo = 0 (i.e. A = (k — 6)/k) it 
is deff = 3 — 8/k. In the physically interesting parameter range k < 8, 
the former is > 2 and the latter is < 2. Recall now that a Bessel process 
is recurrent (not recurrent) if its effective dimension is less (greater) than 
2. Thus, the driving processes X^ ' hit each other almost surely in the case 
hoo = 0 and they don’t hit (a.s.) in the case hoo = h 2 {K). Since the hitting of 
driving processes means hitting of the SLE traces, this teaches us that case 
hoo = 0 describes a single curve joining Xi and X 2 while case hoo = ^ 2 (k) 
describes two curves converging toward infinity. 

Notice that previous results are independent of ai and 02 , provided their 
sum does not vanishes. We also observe that setting ai = 1 and 02 = 0 
(or vice versa) one recovers an SLE(«;;kA). Recall that if hoo = 0 then 
p = kA = k — 6 corresponds to an ordinary chordal SLE from Xi to X 2 . 
Our double SLEs with hoo = 0 corresponds to one chordal SLE seen from 
both ends and the fact that the tips of the traces hit is natural. The other 
case, hoo = ^ 2 (^) corresponds to /? = kA = 2 and since the driving processes 
can not hit, the process can be defined for all t > 0. Assuming that f^(ai + 
a 2 )dt = 00 , the capacity of the hulls grow indefinitely and (at least one of) 
the SLE traces go to infinity. 

The two possible geometries are illustrated in figure 0] 
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8.1.2 A mixed case for 2SLE 


Because of its simplicity, we use double SLE as a testing ground for mixed 
correlation functions. So we consider the sum 


Z — \Zq + 11Z2 


with both A and /x positive. As already mentioned, the interpretation of Z as 
the continuum limit of partition functions of lattice models is unclear since 
Zq and Z 2 do not scale the same way. We nevertheless study it to illustrate 
ways of computing (arch or geometry) probabilities. As one may expect, 
we no longer have an almost sure global geometry but rather nontrivial 
probabilities for the two geometries; either no curve at infinity or two curves 
converging there. 

Let r = inf{t > 0 : be the stopping time which indicates 

the hitting of the driving processes - and thus of the two curves. We can 
define the driving processes as solutions of the 2SLE system on the (random) 
time interval t G [0, r). At the stopping time we define friz) = lim^i,- fsiz) 
for such z G El that the limit exists and stays in the half plane H. The hull 
Kr is defined as the set where the limit doesn’t exist or hits dEI. 

The question of geometry is answered by the knowledge of whether the 
two traces hit, that is whether r < 00 or not. Thus we again consider the 
difference L) = — X^‘^\ whose Ito derivative is now ; 

dEi = ^dBt + |.(ai + a2)dt + ~ ^ (ai + a2)dt 

with Bt = — y/a^Bj.^'^ is a Brownian motion, so that after a time 

change, ds = (ai + a 2 )dt, the result doesn’t depend on oi or 02 . The last 
drift term comes from the derivative of logZ. 

One might for example try to find the distribution of r by its Laplace 
transform KYQ=y[e~^'^] = fpiv)- By Markov property, 


is a closed martingale on t G [0, r) so requiring its Ito drift to vanish leads 
to the differential equation 


(-^ ^ 

V 01 + 02 


{k 


- 6)A + 2/iy(8-0/« 
A + /xy(8-'")/'" 


)9y+'^d‘^)ff3iy) 


= 0 
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The result depends only on /3/(oi + 02 ). We conclude that the distribution 
of (oi + a 2 )r, the capacity of the final hull Kr, is independent of the speeds 
of growth oi and 02 - Also the result depends on A and /r only through ^/\. 

In particular we want to take /3 | 0 to compute the probability that the 
traces hit. Constant functions solve the differential equation but another 
linearly independent solution has the correct boundary values /o( 0 ) = 1 and 
/o(oo) = 0 , namely 


Pyo=y[r < 00 ] = hmEyo=y[e 


A 

A + 


As expected on general ground, this is the fraction of the two partition 
functions AZq and Z = XZq + ^Z 2 . 


8.2 Triple and/or quadruple SLEs 

We will give a few more of examples of multiple SLEs. Certain triple and 
quadruple SLEs are the scaling limits of interfaces in percolation and Ising 
model with rather natural boundary conditions. These models will be con¬ 
sidered in section 10 Here we study triple and quadruple SLEs for their 
own sake. We restrict ourselves to k < 8 . 


8.2.1 3SLE (pure) configurations 

Partition functions with n = 3 have only two possible scaling behaviors 
depending whether the weight h^o of the field at infinity equals either to 
hsiK) = ^ or to /ii(k) = This follows from CET fusion rules. For 

reasons already explained we shall not mixed them. 

The case hoo = h'i{K) is the simplest. There is only one possible partition 
function with this scaling, namely 

[(X 2 - Xi)(X3 - Xi)(X3 - X 2 )f^ 

It is expected to correspond to configurations with three curves starting at 
initial positions Xi, X 2 and X 3 and converging toward infinity. 

The case hoc = /ii(k) is more interesting since the space of such partition 
functions is of dimension two and coincides with the space of conformal block 
with 4 insertions of boundary operators 'll), with one localized at X 4 = 00 : 

(V’(X4)V’(X3)V’(X2)V’(Xi)) 

We assume the points to be ordered Xi < X 2 < X 3 < X 4 . The associated 
process should describe a family of two curves joining any pair of adjacent 
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Figure 5 : For hoo = hi{K) the curves of 3 SLE join either [X1X2] and [X3X4] 
(on the left) or [X4X1] and [X2X3] (on the right). 

points without crossing. There are thus two possible topologically distinct 
geometries: either the curves join the pairs [X1X2] and [X3X4] or they join 
[X4X1] and [X2X3], see figure [SJ As expected, the number of topologically 
distinct configuration equals that of conformal blocks, namely two. Notice 
that the last process is the same as a 4 SLE but with the speed 04 vanishing, 
see figure ini 

By conformal invariance we may normalize the points so that Ai = 0 , 
X2 = X, X3 = 1 and A4 = 00 with 0 < a: < 1 . We have two distinct 
topological configurations and we thus have to identify the two correspond¬ 
ing pure partition functions. This is will be done by specifying the way the 
partition functions behave when points are fused together. By construction 
these partition functions may be written as correlation functions 

Z{x) = (V^(oo)V’(l)V'(x)V^( 0 )) 

so that their behavior when points are fused are governed by CFT fusion 

ht—Q 2 

rules. As a consequence, Z{x) behave either as x~^ or as as x —> 0 . 



Figure 6: Arch configurations for four SLE processes in an arbitrary domain. 
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We select the pure partition functions Zi and Zn by demanding that; 

Zi{x) = x~^ x[l+ •••], 

= (1 — X [const. + • • 

and Zn{x) = Zi{l — x) so that 

2 

Zii(x) = X'^ X [const.+ ••• ], 

. . ft—6 _ - 

= (1 - x) « X [IH-], 

Zj will turn out to be the pure partition function for configurations in which 
the curves join the pairs [Ox] and [loo] while Zj/ will turn out to correspond 
to the configurations [xl] and [ooO]. The rationale behind these conditions 
consists in imposing that the pure partition function possesses the leading 
singularity, with exponent (6 —k)/k, when x is approaching the point allowed 
by the configuration but has subleading singularity, with exponent 2 /k, when 
X is approaching the point forbidden by the configuration. 

This set of conditions uniquely determines the functions Zj and Zj/. 
These follows from CFT rules but may also be checked by explicitly solv¬ 
ing the differential equation that these functions satisfy. Writing Z(x) = 
x^/^(l — x'f‘1'^ G{x) yields, 

K^x{l — x)G''{x) + 8«;(1 — 2x)G'(x) — 4(12 — k)G{x) = 0 
so that G{x) is an hypergeometric function and 

Zjj{x) = const.x^/'^(l — x)^/^ T(—,- -] —Jx) 

K K K 

with the constant chosen to normalize Zj as above. Using this explicit 
formula one may verify that Zj{x) is effectively a positive number for any 
X G [0; 1] so it has all expected properties to be a pure partition function. 
For K = 4, Zj{x) = y/{l — x)/x and for k = 2, Zi{x) = (1 — x^)/x^. 

8.2.2 Arch probabilities 

Let us now compute the probabilities for having one of the two topologically 
distinct configurations; either (/) with curves joining either [Ox] and [loo] 
or {II) with curves joining [xl] and [ooO] as we just discussed. We shall pro¬ 
ceed blindly, but the reader should beware that there are subtleties involved. 

ii) 

What is computed is the probability for certain ’s to hit each other. 


as X ^ 0 (9) 

■ • ], as X ^ 1 

as X ^ 0 
as X ^ 1 
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What happens at the level of hulls and how the process should be prop¬ 
erly continued is not investigated, but is expected to yield the announced 
probability for arch configuration. 

We consider a generic partition function Z which is the sum of the pure 
partition functions Zj and Zu'. 

Z(x) =pjZj(x) +pj/Zjf(x) 

with pj and pu positive. To specify the 3SLE (or 4SLE) process we need 
the partition function Z(Xi, X 2 , X 3 , X 4 ) which is recover from Z{x) by con¬ 
formal transformation : 

Z{Xi,X 2 ,X^,Xi) = [(X 4 -X 2 )(X 3 -Xi)]^ Z{X) 
with X the harmonic ratio of the four points Xi, X 2 , X 3 and X 4 : 

“ VXi -X 3 AX 4 -X 2 / 

Let Mi{x) and Mjj{x) = 1 — Mi{x) be defined by 

Mj{x) = piZj{x)/Z{x) , Mii{x) = piiZii{x)/Z{x) 

By construction the processes t M[{Xt) and t Mjj{Xt), with Xt the 
harmonic ratio of the four moving points, are local martingales. Since both 
Zj and Zjj are positive, Mj are Mjj are bounded local martingales and thus 
are martingales. 

Let r be the stopping time given by the first instant at which a pair 
of points x'f^ coincide. Then, in configuration (I) we have lim ^Xt = 0 
while limt/'rXt = 1 in configuration (II). Since, for k < 8 , Mi{x) is such 
that Mi{x) = 1 but lima,^! Mi{x) = 0, we obtain that Mj evaluated 

at the stopping time r is the characteristic function for events with the 
topological configuration (/), ie: 

limAf/(^j) lconfig.(/) 

liip. Mjj(^Xt) lconfig.(//) 

t/r 

Since Mj and Mjj are martingales, we get the probability of occurrence of 
configurations of topological type (I): 

P[config.(/)] = Mi{Xt=o) = — ^ (10) 

piZi(x) +piiZn{x) 

and similarly for the probabilities of having configuration {II). As expected 
they are ratios of partition functions. 
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8.3 Applications to percolation and Ising model 

We are now ready to give an application of triple (or quadruple) SLE to 
percolation and Ising model. Exploration processes in critical percolation 
are described by SLEs with k = 6, as proved in m- Interfaces of spin 
clusters in critical Ising model is believed to correspond to k = 3 while 
interfaces of Fortuin-Kasteleyn clusters ~ which occur in a high temperature 
expansion of the Ising partition function - are expected to correspond to 
the dual value k = 16/3. 

What we have in mind are these statistical models, defined on the upper 
half plane, with boundary condition changing operators at the four points 
0, X, 1 and oo. They change the boundary condition from open to closed 
(or vice versa) in percolation (k = 6) and from plus to minus (or vice versa) 
for Ising model (k = 3). 

To apply previous results on 4SLE processes to these situations, we have 
to specify the partition functions Z{x), or equivalently, we have to specify the 
value of Pi and pij. This is done by noticing that these models are left-right 
symmetric so that for x = 1/2 there is equal probability to find configuration 
(/) or {II). Since Z/(l/2) = Z//(l/2), we have pi = pn = 1, so that the 
total partition function is Z{x) = Zi{x) + Zii{x) and the probability of 
occurrence of configuration (I) for any 0 < x < 1 is now: 

P[config.(/)] = , . , Zii{x) = Zi{l - x) 

Zi{x) + Zu{x) 

- Percolation corresponds to k = 6. The boundary changing operator 
ijj has dimension 0. The pure partition function Zj has a simple integral 
representation: 

= £d,»-V3(l_,)-2/3. 

By construction Zii{x) = Zi{l — x) also possesses a simple integral repre¬ 
sentation but, most importantly, it is such that the total partition function 
is constant, Z(x) = Zi(x) + Zii(x) = 1, as expected for percolation. As a 
consequence we hnd: 

P[config.(/)]perco = ^( 1 / 3)2 j 
This is nothing but Cardy percolation crossing formula. 
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~ Ising spin clusters correspond to k = 3. The boundary changing oper¬ 
ator il) has dimension 1/2 and may thus be identified with a fermion on the 
boundary. However the pure partition functions do not correspond to the 
free fermion conformal block. By solving the differential equation with the 
appropriate boundary condition we get: 


^7(3:) spin Ising 


1 — X + 

= const- - -^ 

X{1 — X) 



{y{l - y)f/^ 

(1 -y + 2/^)^ 


The total partition function Zi{x) + Zi{l — x) is proportional to 
which is the free fermion result. 

Hence, the Ising configuration probabilities are : 


]P[config. (/)]spjn Ising — 



(1/(1 - 

(1 - ?/ + y2)2 


/■\ (1/(1 - 

/o \l-y + y^Y 


This is nothing but a new - and previously unknown - Ising crossing formula. 

- FK Ising clusters correspond to k, = 16/3. The operator 1 / has then 
dimension 1/16. The pure partition function are given by: 


Zi{x) FK Ising 


(1 — 

xV8(l + ^)l/2 


and the crossing probabilities by: 


P[config.(/)]FK Ising 


Y^(l — x) -|- (1 — 

\/ X + x^/^ + Y^(l — x) + (1 — x)3/2 


The other critical random cluster (or Potts) models with 0 < Q < 4 
have Q = 4cos^ (^),4 <k< 8 and it is straightforward to obtain explicit 
crossing formulas involving only hyper geometric functions. 


8.4 nSLEs and beyond 

We now comment on how to compute multiple arch probabilities for general 
nSLEs. This section only aims at giving some hints on how to generalize 
previous computations. So it shall be sketchy. It is clear that the key point 
is to identify the pure partition functions - once this is done the rest is 
routine. As exemplified above by eq. this is linked to CFT fusions. The 
rules there were that, for a given arch system, fusing two points linked by an 
arch produces the dominant singularity which means that the two boundary 


36 













Figure 7: Illustration of the fusion rules corresponding to arch configura¬ 
tions. 

operators are fused on the identity operator, whereas fusing two points not 
linked by an arch produces the subleading singularity which means the fusion 
of the two boundary fields on the identity should vanish. In general there 
could be a whole hierarchy of arches with arches in the interior of others, i.e. 
with a family of self-surrounding arches, the next encircling the previous. 
So we are lead to propose the following rules. 

For a given arch configuration: 

— The most interior pair of adjacent pair of points, say Xj, Xj+i in a family 
of self-surrounding arches fused into the identity operator, so that the pure 
partition function evaluated at Xi ~ Xj+i should be proportional to (Xj+i — 

K. — 6 

Xi)~^ times the pure partition function associated to the arch system with 
the interior arch [XjXj+i] removed. Symbolically : 

Zpurei' ■ ■ ) Xi — Xj_|_i, •••);:; const. (Xj_|_i — Xi) « X .^pure\[XiXi+i](' ■ ■ ) ■ ■ ■ ) 
for Xi and Xj+i linked by an arch. 

— The fusion on the identity of any pair of adjacent points not linked by 
an arch should vanish, so that the fusion of this pair of points produces the 
subleading singularity. Symbolically : 

2 

•^pure(' ■ ■ ) Xi — Xj_|_i, • • • j ~ const. (Xj_|_i Xi) K • 

for Xi and Xj_|_i for not linked by an arch. 

We do not have a complete proof that these rules fully determine the 
pure partition functions but we checked it on a few cases, see figure (Tj 

Here are a few samples. We shall give the relation between the pure parti¬ 
tion and the CFT conformal blocks indexed by the corresponding Bratelli di¬ 
agram. For n = 4, we may have the following arch systems [XiX 2 ][X 3 X 4 ] or 
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[Xi[X 2 X 3 ]X 4 ]. (A given geometrical configuration may correspond to differ¬ 
ent arch systems depending at which location we open the closed boundary. 
But they are all equivalent to these two up to an order preserving relabeling 
of the points. For instance [X 4 A 1 ][X 2 A' 3 ] is equivalent to [Xi[X 2 X'i]X^.) 
Applying the previous rules we get: 

^[Xi[X2X3]X4] = 

^[XiX2][X3X4] = {[ho]i'{Xl)[h^]i'{X2)[ho]i^{X^)[hi]'4’{Xi)[ho\) 

where the indices /im, m = 0,1, • • • refer to the corresponding points in the 
Bratelli diagram, i.e. to the weights hm{n) of the intermediate Virasoro 
modules. The coefficient u is fully determined, in terms of CFT fusion 
coefficients, by demanding that the fusion of X 2 and X 3 on the identity 
vanishes. 

One may go on and solve for the pure partition functions in few other 
cases. A particularly simple example with n = 6 is given by : 

Z[Xi[X2[X3Xi]X5]Xe] = 

{ [ho] V’ ) [hi] V' (-^^2) [fe 2 ] V' (^3) [h3] V’ {Xi ) [h2] V’ (A' 5 ) [hi] V' (-^^6 ) [ho]) 

As can be seen on these examples, there is no simple relation between arch 
systems and Bratelli diagrams and the change of basis for one to the other 
is quite involved. The only simple rule we find is that the pure partition 
function for a unique family of self-surrounding arches is a pure conformal 
block corresponding to a unique Bratelli diagram. 
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